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ABSTRACT 

We present optical and near-infrared light curves and optical spectra of SN 2013dx, associated with 
the nearby (redshift 0.145) gamma-ray burst GRB 130702A. The prompt isotropic gamma-ray energy 
released from GRB 130702A is measured to be = 6.4^1;g x 10®° erg (IkeV to lOMeV in the 

rest frame), placing it intermediate between low-luminosity GRBs like GRB980425/SN 1998bw and 
the broader cosmological population. We compare the observed g'r'i'z' light curves of SN 2013dx 
to a SN 1998bw template, finding that SN 2013dx evolves ~ 20% faster (steeper rise time), with 
a comparable peak luminosity. Spectroscopically, SN 2013dx resembles other broad-lined Type Ic 
supernovae, both associated with (SN 2006aj and SN 1998bw) and lacking (SN 1997ef, SN 20071, and 
SN 2010ah) gamma-ray emission, with photospheric velocities around peak of ~ 21,000kms“®. We 
construct a quasi-bolometric [g'r'i'z'yJ) light curve for SN 2013dx, only the hfth GRB-associated SN 
with extensive near-infrared coverage and the third with a bolometric light curve extending beyond 
At > 40 d. Together with the measured photospheric velocity, we derive basic explosion parameters 
using simple analytic models. We infer a ®°Ni mass of Mni = 0.37 ± 0.01 Mq, an ejecta mass of 
Mej = 3.1 ±0.1 Mq, and a kinetic energy of Ek = (8.2±0.43) x 10®® erg (statistical uncertainties only), 
consistent with previous GRB-associated supernovae. When considering the ensemble population of 
GRB-associated supernovae, we find no correlation between the mass of synthesized ®®Ni and high- 
energy properties, despite clear predictions from numerical simulations that Mni should correlate with 
the degree of asymmetry. On the other hand, Mni clearly correlates with the kinetic energy of the 
supernova ejecta across a wide range of core-collapse events. 

Subject headings: gamma-ray burst: individual (GRB 130702A) — supernovae: individual (SN 2013dx) 
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1. INTRODUCTION 

The evidence for the association between long-duration 
gamma-ray bursts (GRBs) and the death of massive 
stars has been steadily growin g o ver the last two decades 
(see Woosley & Bloom) 2006 and Hjorth & Bloom 2012 
for reviews]. 'i'he hrst direct evidence of this link was 
a spatially and temporally coincident supernova (SN), 
SN 1998bw, with GRB 980425 at redshift z = 0.0085 


( Gala ma et al.|1998||Iwamoto et al.|1998[|Kulkarni et al. 

iJ M Aiu ,UL _ II _ L^_ _ L ^ __ „l _ 


yysp. Since dlN lyOsbw, there have been a number of 
spectroscopically confirmed supernovae (SNe) associated 
with GRBs (Table [^. 

While most, if not all low-z long-duration GRBs ap¬ 
pear to be accompanied by SNe (the exceptions being 
GRB 060614 and GRB 060505; IFynbo et al.||2006| |GaH 
Yam et al. 2006 Della Valle et al. 2006), only a small 
traction of core-collapse e xplosions are capable of gen¬ 
erating relativistic ejecta (iBerger e t al.||2003[ [Soderberg 
et al. 2010 Bietenholz et al. 2014). Even when limited 
to the specihc subtype of SiNe associated with GRBs, 
the broad-lined Type Ic SNe, those with and without 
relativistic ejecta appear to be indistinguish able based 


on their light curves (e.g., Drout et al. 2011). However, 


spectra of the host galaxies reveal that GRB-SNe pre¬ 
fer more metal-poor enviro nments than Type Ic-BL SNe 


without asyciate d GRBs (jModjaz et al. 
&: Fruchter||2013 1. 


2008 Graham 


I'urthermore, within the GRB population, there is a 
considerable diversity in the observed prompt gamma- 
ray energies spanning six orders of magnitude, = 

1048_io54 g]-g_ suggested that low-luminosity 

GRBs {Ej^iso ^ 10^® erg) have “failed” jets that can¬ 
not pierce their stellar envelope and instead dissipate 
energy into the star to create relativistic shock breakout 


(Bromberg et al. 2011 


Margutti et al. |m4l|Nakar|m5| ). 
different appearance at high ener- 


But despite their very different appearance at hig. 
gies, as of yet there is no clear distinction between SNe 
associated with low-luminosity GRBs (e.g., SN lOOShw) 
and the larger (observed) cosmological population (e.g., 
SN 2003dh, SN 2013cq). 

With still only a handful of well-observed examples, 
each new nearby GRB affords a unique opportunity to 
understand the central engine powering these outflows. 
In particular, we can probe the progenitor from two dif¬ 
ferent angles by studying the SN simultaneously with 
the GRB. Here we present observations of SN 2013dx as¬ 
sociated with GRB130702A. At z = 0.145, SN 2013dx 
is sufficiently nearby to enable a detailed photometric 
and spectroscopic study of the SN evolution. Further¬ 
more, with Ej^iso = 6.4 X 10®° erg (for IkeV to lOMeV 
in the rest frame), the prompt-emission properties place 
GRB 130702A between most low-luminosity GRBs and 
the more energetic cosmological population. 

Throughout this paper we use the convention F^{t) oc 
photon index T = ,0-1-1. We assume a AGDM 
model with Hn = 60.6 km s~® Mp c~^, 0^ = 0.286, and 
Ha = 0.714 (Be nnett et al. 20 14|). A ll photometry is in 
the AB system (|Uke &: Gunn||iy83|), and quoted uncer¬ 
tainties are Icr (687o) conhdence intervals unless other¬ 
wise noted. Dates and times are UT in all cases. 


GRB 130702A was detected by the Gamma-Ray Burst 
Monitor (GBM; Meegan et al.|2009 1 o n the Fermi satel- 
lite at 00:05:23.070 on 2013 July 2 (|Collazzi & Gon- 
naughton |2013 ). As observed by the GBM, the prompt- 
emission duratioij^ was Tgg « 50s (50-300keV). High- 
energy emission was also detected by the Fermi Large 
Area Telescope (LAT; Atwood et al.||200^ [Cheung et al. 


2013 ), as well as by Konus-Wind (Golenetskii et al. 
20131. We adopted the Konus-Wind 20-1200 keV tlu- 
ence of fj = (6.7 ± 0.8) x 10“®ergcm“^ (assuming a 
power-law spectrum with photon index T = 1.87± 0.11). 
Employing the wide-held imaging and rapid transient 


Transient Factory (iPTF; 

Law et al. 

2009 

Rau et al. 

20091, 

Singer et al. 

(20131 

r-f. TTTT - 

^ discovered the optical at- 


as iPTF13bxl, is located at (J2000.0) coordinates a = 
14‘'29“14".78, 6 = -fl5°46'26"4. 

Subsequently, the redshift of GRB 130702A was deter¬ 
mined to be z = 0.145 based on the detection of narrow 
host-galaxy e mission lines ([O III] and Ha) at the after¬ 
glow location (Mulchaey et al. 2013b|a D’Avanzo et al. 


20131 

cate* 


Several other galaxies at or near this redshift lo- 


in the held indicate that GRB 130702A occurred 
in a group or clu ster environment, which is highly un - 
usual for a GRB (Kelly et al. 2013 D’Elia et al. 12015 ). 
The GRB host galaxy may be a metal-poor satellite of 
an adjacent massive spiral (SDSSJ142914.57-1-154619.3), 
which has an offset of only ^ 19 kpc in projecte d dis¬ 
tance and < 60kms“^ in line-of-sight velocity (Kelly 


At z = 0.145, the observed Konus-Wind huence cor¬ 
responds to an isotropic energy release of E^^iso = 
X 10®° erg (IkeV to lOMeV in the rest frame). 
This places GRB 130702A securely between the low- 
luminosity class of events represented by GRB 980425 
/ SN 1998bw and typical cosmologically distant events 
w ith E^^^n > 10®^ erg . 


Butler et al. (2013) reported a battening of the opti- 
cal afterglow 5.26 d after the burst. A spectrum taken 
^ 6 d after the burst showed broad featu r es resembling 
those of SN 1998bw (ISchulze et al.||2013|). |Cenko et al. 
(|2013|) and iD'Elia et al.| (|2013| obtained spectra that 
conhrmed the presence of an emerging SN, dubbed SN 
2013dx, and i dentihed simil a rities with SN 1998bw and 
SN 2006aj. D’Elia et al. (2015) (hereafter D15) re¬ 
ported GRB l30702A/SiN 2()13dx light curve, spectra, 
and SN energetics properties with which we will com¬ 
pare throughout this paper. 

3. OBSERVATIONS AND DATA REDUCTION 
3.1. Optical/Near-Infrared Imaging 

We obtained optical and near-infrared (NIR) images 
of the location of GRB 130702A / SN 2013dx with the 
robot ic Palomar 60 inch telescope (P60; i Cenko et al. 
' 20061 ), the 2 m Liverpool Telescope (LT; |Steele et al. 
20041, the Reionization and Transients Infrared/ (Jptical 
camera on the 1.5 m Harold L. Johnson Telescope 


(RATIR; Butler et al. 2012 Watson et al. 2012), the 
Large Monolithic Imager on the 4.2 m Discovery Channel 


2. GRB130702A/SN 2013DX 


Tgo is defined as the time over which a burst emits from 5% 
of its total measured counts to 95%. 
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TABLE 1 

GRB-SN References 


GRB-SN 


References 


GRB 980425/SN 1998bw 
GRB 030329/SN 2003dh 
GRB 031203/SN 20031w 
GRB 060218/SN 2006aj 

GRB 091127/SN 2009nz 
GRB 100316D/SN 2010bh 

GRB 120422A/SN 2012bz 
GRB 130427A/SN 2013cq 
GRB 130702A/SN 2013dx 
GRB 140606B/iPTF14bfu 


Galama et al. 


Stanek et al.| 


UD^t 


iP 


19981: llwamoto et al. 


Malesam et al'!TJ:^uu4TT T Gal-Yam et al. 
Oampana et al'l i ^UUb|i ; ^Iodiaz et aL| 2uub ^ 
soiierman et ai.l | 2 UUt)ll;fi-'ian et aij dzCroHlrTn^ 


latheson et al. 


1 1998 ; [Kulkarni et al.| l |1998[ | 


2UU3 I 


Cobb et al, 


Starling et aU "f 


uano et~^ 


Ivlelandri 


et ai.j | 2 uun||;^ian et ^ jTl u^;| 

B Z Ululn' B erger et al.|l|2 3 iili 
T 1f2rr i; Tl; iBufano et aP l|2Ul2 l; 


Thomsen et al. 


Vlirabal et al. 


CocevsKi et ai. 


unli ; IFerrero et al.M2006ll 

UD^ '-' 


j"^dSchulze et a^ l l|2Ui4|| 
i jzuij i; [Levan et al.| (['2U14[ |; [Melan 


ChornocK et ailJ^uiu 


; [OlI 


livares E. et al 




Xu et al. 

mis wor 


Cano et al. 


rmia et all 

|2Ulb|| - - 


[Melandri et al.| (|2014^; [Perley et al.|(|201^ 


Note. — We do not include the recent detection of GRB 111209A/SN 2011kl associated with a superluminous SN 
( [Greiner et al.|2015[ | because it is believed to be powered by a magnetar and not solely powered by ^®Ni. 


TABLE 2 

Log of Spectroscopic Observations 


At 

(d) 

Instrument 

Exposure 

(s) 

Wavelength Range 
(nm) 

Slit 

(") 

Grating/Grism 

Airmass 

1.17 

DBSP 

1800 

340-1000 

1.0 

1200/5000 -1- 1200/7100 

1.06 

3.25 

DBSP 

1800 

340-1000 

1.0 

600/4000 -1- 316/7500 

1.36 

6.22 

DBSP 

1800 

340-890 

2.0 

600/4000 -1- 600/10000 

1.25 

9.33 

DEIMOS 

600 

450-950 

1.0 

600/7500 

1.19 

11.34 

DEIMOS 

600 

490-1010 

1.0 

600/7500 

1.32 

14.21 

DBSP 

1200 

350-1000 

1.5 

600/4000 -1- 316/7500 

1.27 

31.28 

DEIMOS 

900 

450-950 

1.0 

600/7500 

1.30 

33.27 

LRIS 

1200 

330-1020 

1.0 

400/3400 -1- 400/8500 

1.30 

330.39 

LRIS 

1460 

330-1020 

1.0 

400/3400 -1- 400/8500 

1.04 


Note. — At is the time from Fermi trigger in observer frame. 


Telescope (LMI/DC T), and the Low Resolution Imaging 
Spectrometer (LRIS; |Oke et al. 11995 ) on the 10 m Keck-I 
telescope. Additionally, we included early-time Palomar 
48-inch r' obs ervations (O.U d and 0.21 d after the GRB 
trigger) from Singer et al. (2013a). The reduction pro¬ 
cedures for each individual facility are described below, 
while the resulting photometry is presented in Table 
and plotted in Figure Photometry from different tele¬ 
scopes is calibrated to same stars for uniform calibration. 
After rem oving extinction, afterglow, and host galaxy 
(see § 4.1), the cross-calibration errors are ~ 0.03-0.05 
mag (approximately 3-5% in flux). 


3.1.1. P60 Photometry 

P60 observed the location of GRB 130702A in the g', 
r', i', and z' filters beginning 1.17 d after the Fermi GBM 
trigger. Basic GGD reductions are provided in real time 
by a custom IRAlEj/PyRA Ip] pipeline. A t later times 
images were stacked with SWs 


(At> 3 d), images ’ 


SWarp (Bertin et al. 


2002) on a nightly basis to increase the signal-to-noise ra¬ 
tio (SNR). We performed aperture photometry at the af¬ 
terglow location, calibrating with respect to near by point 
sources from the Sloan Digital Sky Survey (SDSS;|Aihara 

eFaI][MT| ). 


IRAF is distributed by the National Optical Astronomy Ob¬ 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA), Inc. under cooperative agree¬ 
ment with the National Science Foundation (NSF). 

See http://www.stsci.edu/institute/software_hardware/ 

pyraf 


3.1.2. LT Photometry 

LT began observing the location of GRB 130702A with 
the 10:0 GGD camera 9.87 d after the GBM trigger. Ob¬ 
servations were obtained in the g', r', i', and z' filters. 
Standard reduction techniques were applied to detrend 
the data, and photometry was performed in the same 
manner as for the P60 images (including the same SDSS 
reference stars for photometric calibration). 


3.1.3. RATIR Photometry 

RATIR obtained simultaneous multi-color {r'i'z'yJFl) 
imaging of the location of GRB 130702A beginning 2.16 d 
after the GBM trigger. The RATIR data were reduced 
using an automatic python pipeline with bias subtrac¬ 
tion and twilight-sky flat fielding. Given the lack of a 
cold shutter in RATlR’s design, IR dark frames were not 
available. Laboratory testing, however, confirmed that 


the dark cu rrent is negligible in both IR detectors (Fox 
et al.|[2M2). Astrometric solutions were calculated from 


astrometry.net (Lang et al. 2010) and the individual 


frames are stacked using SWarp. 

We performed aperture photomet ry on the resulting 


stack ed images using Sextractor (Bertin & Arnouts 
1996) with an inclusion radius determined from the me¬ 


dian full width at half-maximum intensity (FWHM) of 
the images. The resulting instrumental magnit udes were 
compared t o SDSS in the optical and 2MASS (|Skrutskie| 
et al. 2006) in the NIR to calculate zeropoints. I'or the 
y band, we created a spectral energy distribution (SED) 
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m 



Time since trigger (days) 


Fig. 1.— Optical and NIR light curves of GRB 130702A, corrected for Galactic extinction and host extinction (assuming A^/host = 0.10 
mag). Note that the host-galaxy contribution is not removed. Diamonds are P48/P60 data, circles are RATIR data, and stars are Li verpool 
data: error b ars are overplotted on all datapoints. Solid lines indicate power-law fits with a\ = 0.57 and tf, = 1.17d taken from |Singer| 
|et al.| (|2013a|l and 02 = 1.25 from the XRT power-law decay index. Dashed lines mark the extrapolated power law where we assume 02 
tor all times beyond tf,. 


from the combination of optical and NIR catalog sources 
and interpolated to the appropriate wavelength. To place 
all photometry on t he AB system, we used t he J- and 
iJ-band offsets from Blanton & Roweis (2007). 


3 . 1 . 4 . Keck/LRIS Photometry 

The location of GRB 130702A was observed with 
Keck/LRIS on 2014 May 28 (At = 330 d) in the u', g', 
and i?-band filters. The resulting images were reduced 
using the LPipe packag^^ Because the host galaxy 
was clearly resolved in some of the better-seeing images 
(FWHM = (y.'75), we adopted an aperture radius of I"5 
to incorporate all of the flux from the visible extent of 
the galaxy (Figure]^. Photometric calibration was per¬ 
formed relative to point sources from SDSS. 


3 . 1 . 5 . Keck/MOSFIRE Photometry 


We imaged the location of GRB 130702A with the 
Multi-Object 
(MOSFIRE; 


McLean et al. 

2012) on the 10 m Keck I 

2014 June 16. 

mages were obtained in the 


J and Ks filters and reduced using custom IDL scripts. 
We perform ed ape rture photometry using a 1"5 inclusion 
radius (see § 3.1.4), with photometric calibration relative 
to 2MASS. 


3 . 1 . 6 . LMI/DCT Photometry 

The location of GRB 130702A was observed with 
LMI/DGT on 2015 March 27 (At = 633 d) in the 5 ', r', i', 
and z' filters. The resulting images were detrended with 
a custom IRAF pipeline. Indi vidual frame s were astro- 
metrically aligned with Scamp ( Bertin|2006 ) and coadded 
using SWarp. Photo metry was calculated in the manner 
described in § |3.1.3[ 


See http://www.astro.caltech.edu/~dperley/programs/ 
lpipe.html tor details. 


3 . 2 . Optical Spectroscopy 
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Fig. 2.— Keck/LRIS g'-band image of the host galaxy of GRB 130702A, obtained at At = 330.48 d after the GBM trigger (i.e., when 
the afterglow and SN emission had faded away). The location of the transient is displayed in the inset with the white cross. The dwarf 
host is clearly elongated in the N-S direction, with the bulk of the star formation (as evidenced by the transient location and the nebular 
emission lines) apparent in the northern component. The image is oriented with N up and E to the left. 


We obtained a series of optical spectra of 
GRB 130702A, beginning at At = 1.17d after the 


the XRT Light Curve Repositorjj^ The time-averaged 
spectrum was well described (W-stat = 299.10 for 374 de- 


(DBSP; Oke & Gunn 

1982|) on the 5 m Palomar Hale with no evidence for Ah in 

excess of the Galactic value 

Kalberla et al.||2005 

I. 

telescope, Keck/LKib on Keck-i, and the Dhiep imaging (IVh Gai = x cm 


Mult i-Ob.iect Spectrograph (DEIMOS; Faber et al. 


20031 on the 10 m Keck-II telescope. An observing log 


IS presented in Table All spectra were obtained with 
the slit oriented at the parallactic angle to mini mize 


differential lo sses due to atmospheric dispersion (Fil- 
|ippenko||1982[ though note also that LRIS employs an 
Atmospheric Dispersion Corrector to further mitigate 
against differential slit losses). The resulting reduced 
one-dimensional spectra are displayed in Figure 
All spectra were r educed using standard routines and 
optimally extracted (|Horne|1986|i wit hin the IRAF envi¬ 
ronment (see, e.g., Cenko el al.||2008 for details). A dis¬ 
persion solution was computed using calibration spectra 
of comparison lamps, and then adjusted for each indi¬ 
vidual exposure using night-sky emission lines. For the 
LRIS and DFIMOS spectra, sky background e mission 
was su btracted using the algorithm described by |Kelson| 
(2003). Telluric atmospheric absorption features were 


removed using t he continuum from spectrophotometric 


standard stars (Wade & Horne 1988 Matheson et al. 


2000). Finally, a sensitivity function was applied using 


observations of spectrophotometric standards at a com¬ 
parable airmass. We caution that the final Keck/LRIS 
spectrum obtained at At = 330.39 suffered from a failure 
of the blue shutter, which may impact the relative-flux 
calibration. 

Upon publication of this manuscript, all one¬ 
dimensional spectra will be made publicly available via 
the Weizma nn Interactive Supernov a data RFPository 
(WISeRFP; |Yaron fc GaLY^[M^ . 


GRB 130702A with the Chandra X-ray Observatory on 
2013 September 5, corresponding to At = 65.2 d since 
trigger (PI R. Margutti). The data were reduced with the 
CIAO software package (version 4.6) and corresponding 
calibration files. Standard ACIS data filtering has been 
applied. In 14.9 ks of observations we find clear evidence 
for X-ray emission at the location of GRB 130702A, with 
significance > 50(7. The spectrum was well modeled by 
an absorbed power law with T = 1.66 ± 0.15, consis¬ 
tent with the Swift-XKT time-averaged spectrum. We 
found no evidence for an intrinsic absorption component, 
with a 3(7 limit of fVn.host < 1-5 x 10 ^^cm“^. Adopt¬ 
ing these spectral parameters, the unabsorbed flux is 
Fx = (1.20 ±0.08) X 10-13 erg s-icm-2 (0.3-10 keV). 

A second epoch of Chandra observations was obtained 
on 2013 December 6 (At = 157.5 d) with an exposure 
time of 34.6 ks. GRB 130702A was clearly detected with 
significance > 40(7, which allows us to constrain the spec¬ 
tral evolution (or lack thereof) of GRB 130702A at very 
late times. Our spectral analysis reveals no evidence for 
spectral evolution. The best-fitting power-law index was 
P = 1.85 ± 0.16, with Ah, host < 1-6 x lO^i cm”^ at 3(7 
confidence level. The corresponding unabsorbed flux is 
Fx = (4.5 ± 0.45) X lO”!^ ergs-icm”^ (0.3-10 keV). 

The X-ray light curve of GRB 130702A, comprising 
Swift-XBT and Chandra observations, is presented in 
Figure]^ 

4. LIGHT-GURVE ANALYSIS 
4.1. Isolating the Supernova Component 


3.3. X-Ray Observations 



Emission from the location of GRB 130702A results 
from three distinct components: the GRB afterglow, the 
associated SN, and the underlying host galaxy. Here we 
try to isolate the emission resulting from the associated 
SN, including a proper accounting for line-of-sight ex¬ 
tinction, in order to study the properties of SN 2013dx. 

See http : //vww. swift. ac. u k/xrt curves and the associated 
description in [Evans et al.| ||2UU9|. 



















































3000 4000 5000 6000 7000 8000 9000 
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Fig. 3.— Spectra uncorrected for slit losses, extinction, afterglow, or host-galaxy contamination. Early-time spectra are dominated by 
the afterglow component. The broad features associated with SN 2013dx become visible after about a week. 








7 


D) 

OJ 


X 

3 



Fig. 4.— X-ray light curve of GRB 130702A. 

First, we correct our broadband photometry for ex- 


tinction, both in the Milky Way and 
For the Galactic component, we em 

in the host galaxy, 
ploy the dust-map 

calibration of Schlafly & Finkbeiner ( 

20111 , resulting in 

K(B — V )mw = 0.038 mag, and the 

Milky Way extinc- 


tion law of Cardelli et al. (1989). In order to estimate the 
host extinction, Ay^host, we create an SED at At = 2.25 d 
from linear interpolation. We assume the observed emis¬ 
sion at this stage will be dominated by the (synchrotron) 
afterglow, and thus we fit the SED to a simple powe r-law 
model of the form f^, oc (e.g.. Sari et aL| 1998). We 
incorporate Tly^host as a free parameter, assum ing a Smal l 
Magellanic Cloud (SMC)-like extinction law (Pei 1992). 
We find Ay = 0.13±0.23mag and j3 = 0.52±(J.19 with a 
reduced Xred — 0-83. Adopting Large Magellanic Cloud 
(LMC) and Milky Way (MW) dust extinction laws did 
not alter the derived parameters or fit quality. Eor the 
rest of this paper we assume Ay^host = 0.10 mag. This is 
consistent with other GRB host-extinction values: 50% 
of GRBs have A yhost, < 0-4 mag and 87% of GRBs have 
Ay host < 2 mag (jCovino et al.||2013|). 

Next, we attemipt to remove a ny contribution from the 


afterglow. Singer et al. (2013a I modeled the early-time 


optical emission (At ^ 4 d) as a broken power law with 
an initial decay index of ai = 0.57±0.03 up to the break 
time, tb = 1A7 ± 0.9 d, after which the model followed a 
power-law decay index of a 2 = 1.05 ± 0.03. We repeat 
this analysis with our larger photometric dataset and find 
mostly consistent results. However, even as early as sev¬ 
eral days post-trigger, the observed emission will likely 
have some contribution from the emerging SN (e.g., the 
broad features becoming ^parent in the At = 3.25 d 
DBSP spectrum in Figure!^. Gonsequently, the true af¬ 
terglow decay index may be steeper than indicated here. 

Instead, we consider the decay of the corresponding X- 
ray emission, which is unlikely to be contaminated by the 
SN at At > 1 d. We fit the X-ray light curve to a power- 
law model and find ax = 1-25 ± 0.03. Gombining this 
with the measured X-ray spectral index from the Swift- 
XRT data, I3x = 0.84 ± 0.1 2, we can use standa rd after¬ 
glow closure relations (e.g., |Racusin et al.||2009l and ref¬ 
erences therein) to evaluate where the X-rays fall on the 
broadband synchrotron spectrum. The best fit is found 


for a constant-density circumburst medium with the X- 
rays falling below the synchrotron cooling frequency, t'c. 

As a result, the optical bandpass must fall below t'c 
as well, which is consistent with the measured X-ray to 
optical spectral index of (iox ~ 0.7 at At = 2d. Thus, 
since the optical emission falls on the same segment of the 
SED as the X-rays, it should decay with the same power- 
law index, ao = 1.25. We further assume that both the 
optical afterglow spectral and temporal indices remain 
constant in time over the course of our observations. We 
use this model to calculate the afterglow contribution for 
all our photometric observations. 

We note that this is a significantly shallower decay in¬ 
dex than the ao = 2.2 adopted by D15, who did not 
incorporate multi-wavelength observations into their af¬ 
terglow analysis. A steeper ao may overestimate the 
SN flux at early times, but is negligible at the peak and 
late times when the SN is significantly brighter than the 
afterglow. 

Einally, we must remove the contribution from the un¬ 
derlying host galaxy. In our best-seeing images at late 
times (FWHM « (f.'Tb), the host is clearly resolved, with 
the afterglow/SN location falling on a blue “knot” to the 
north (Eigure0. This location is also responsible for 
the nebular emission lines seen in the final Keck/LRIS 
spectrum (Eigurej^. 

To remove the nost contribution from our photome¬ 
try, we adopt host flux values from our late-time DCT 
{g'r'i'z') and Keck/MOSEIRE (J) imaging and directly 
subtract these from the measured transient fluxes. While 
this does not account for the resolved nature of the host, 
since the typical seeing in our SN data (~ 1'.'5) is com¬ 
parable to the size of the extended emission, this should 
have minimal impact on the resulting photometry. We 
do not have host-galaxy detections in the y and H bands, 
but the host contribution is negligible compared to the 
afterglow and SN components. 

The resulting SN light curves are displayed in Fig¬ 
ure The peak times of the light curves are useful for 
constraining theoretical models, in particular the convo¬ 
lution of total ejected mass, kinetic energy, and opac¬ 
ity of the SN explosion (via the diffusion time). We 
measure the rest-frame peak times by fitting a second- 
order polynomial at 7 < At < 20d and find the fol¬ 
lowing: tp{g') = 11.7 ± 0.3d, tp(r') = 13.2 ± 0.3d, 
tp(j') = 14.7 ± 0.6 d, and tp(z') = 15.1 ± 1.6 d (statis¬ 
tical uncertainties only). These values generally agree 
well with those reported by D15. We do not include 
peak times for yJH because the data are not well sam¬ 
pled close to the peak. 

There is marginal evidence in the i'-band light curve 
(Figure]^ of a decline in flux at early times (At < 3d). 
This is consistent wi th early signatures of shock cool- 
ing (e.g., SN 2006 ai, Campana et al. 2006[ SN 2010bh, 


Gano et al. 2011). However, shock breakout should be 
significantly stronger in bluer bands and we see no indi¬ 
cation of it in either the g' or r' bands. The relatively 
bright optical afterglow of GRB 130702A, compared to 
the optical afterglows of (for example) GRB 060218 and 
GRB 100316D, greatly complicates isolating the SN com¬ 
ponent at early times. Thus, it is difficult to reach firm 
conclusions regarding the presence or absence of a shock- 
breakout signature. 
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Fig. 5.— SN 2013dx spectra with host galaxy and GRB afterglow removed. They are smoothed using a Savitzky-Golay filter with a 30 A 
window. We exclude d th e spectrum at t = 1.17d because it is dominated by the afterglow. Red filled circles mark the position of Si II 
A6355 calculated in § |5.3| used to determine the photospheric velocity of the ejecta. The uncertainties are smaller than the symbol width; 
they are 20 A for t = 9.33 d and 10 A for subsequent epochs. 
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Fig. 6.— Observed SN 2013dx g'r'i'z' data (top left, top right, bottom left, and bottom right, respectively), with synthetic SN 1998bw 
light-curve fit (dashed red li ne) a nd an optimized synthetic SN 1998bw light-cu rve fit scaled by a peak-a r nplitude factor fe and ti me-stretch 
factor s (so lid re d line; see ^4.2|| ov erlaid. SN 1998bw template created from |Clocchiatti et al.| (|2011||, |Galama et al.| ||1998|l, |Sollerman| 
|et al.| l [2002[ |, andjPatat et ai‘] || 1^001 [l. Diamonds are P48/P60 data, circles are KAilK data, and stars are Liverpool data. 


TABLE 3 

SN 1998BW Template Fits 


Filter 

s 

k 

X?ed 

g' 

0.76 ± 0.05 

1.02 ± 0.06 

7.78 

r' 

0.79 ± 0.01 

1.05 ± 0.02 

4.67 

i' 

0.82 ± 0.02 

0.92 ± 0.03 

4.07 

z' 

0.74 ± 0.04 

0.91 ± 0.05 

2.85 


4.2. Comparison with SN 1998bw 

Following past studies of GRB-associated SNe in the 
literature, we next attempt to compare SN 2013dx to the 
well-studied SN 1998bw (associated with GRB 980425). 
We create K-corrected synthetic SN 1998bw light curves 
in the g'r'i'z' filters at the redshift of SN 2013dx, z = 
0.145, using methods described by |Hogg et al. ( 2002^. We 
utilize SN 1998b w photometry anospec t ra from Glocchi- 
atti et al. (2 011|), Galama et al. (19981, Sollerman et al. 
12002]), and IPatat et ai.| (|2001|) I'he K-corrected syn¬ 
thetic SN 19y8bw light curves were also time dilated to 
match the observer frame of SN 2013dx. 

Owing to gaps in the temporal coverage of SN 1998bw 
photometry, especially in the rising phase, we fit the syn¬ 
thetic SN 1998bw light curve in each filter with the em¬ 


pirical functional form from Gano et al. (2011), 


U{t) = A + pt 


/ e(-‘“VF) \ 

yl -I- e(P -*/«)) 


-bt“Rog(t-“3), 


( 1 ) 


allowing A,p, F, R,p,ai,a 2 , and to vary. The result¬ 
ing SN 1998bw synthetic light curves are plotted as 
dashed lines in Figure 

We then assume the hght curves of SN 2013dx in each 
of our four filters can be modeled by simply varying the 
peak amplitude (fc) and stretch factor (s): 

Ll3dx{t) = kUgshw(t/s). (2) 

The resulting fits are plotted in Figure while the mea¬ 
sured stretch and amplitude values are displayed in Ta¬ 
ble We find that SN 2013dx has a peak flux compara¬ 
ble to that of SN 1998bw in all four filters reported here 
(slightly more luminous in the bluer filters, and slightly 
less luminous in the redder filters). With a constant 
stretch value of s « 0.8 in all four filters, the evolution 
of SN 2013dx (in particular the rise time) is noticeably 
faster than that of SN 1998bw. 

However, it is also clear from the fits that SN 1998bw 
is not an ideal match to SN 2013dx, especially in the 
redder i- and z-band filters. Given these important dif¬ 
ferences, we refrain from drawing any physical inferences 
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(e.g., Mni, i?K) from these fits, and instead use the more 
model-indepen dent bolometric light curv e in § [^ 

Cano (2014) and Lyman et al.| (| 2014a l suggeste d that 


GRH-SlNe can be used as standardizable candles. ICanol 


(j2014|) proposed that s and k (measured relative to SN 
lyybbw) are analogous to the absolute peak SN magni¬ 
tude and the amount the light curve fades from maxi¬ 
mum hghCtolSd later (Amis) used in the Phillips rela¬ 
tion ( Phillips|1993 ) for SNe la. We add the uncertainties 
from our s and k ht in quadrature and find that our mea- 
surements of SN 2013dx dev iate by 3.7a from the |Cano 


(2014) hts. The worst fits to 
r' at 3.7a\ however, z' is wit 


Cano (2014) are for g' and 
lin 0.9tT. 4'his further sup¬ 
ports the notion that SN 1998bw is not a perfect match 
for all CRB-associated SNe. 

4.3. Bolometric Light Curve 

We construct the quasi-bolometric light curve of SN 
2013dx using our photometry in the g', r', z', z', y, and 
J filters [H had only upper limits at At > 5d). We in¬ 
clude synthetic photome try for At = 31.28 and 33.27d 
from our spectra (see § 5.1) to supplement photomet¬ 
ric coverage at these epochs. We assume a 10% flux 
error on all synthetic photometry data points. We con¬ 
vert the extinction, host-galaxy, and afterglow-corrected 
magnitudes to monochromatic fluxes and create SEDs 
from linear interpolation of the data for each epoch be¬ 
tween 1 d and 70 d with 0.25 d spacing. Epochs that are 
> 0.5 d from observations are removed to mitigate linear- 
interpolation errors, which only affects observations in 
the y and J bands (see below). We assume that the flux 
is constant across the bandwidth of each filter and use 
trapezoidal integration to calculate the quasi-bolometric 
luminosity. We note that our photometry provides cover¬ 
age over the observer-frame bandpass 0.4-1.35/zm (rest- 
frame 0.35-1.18/rm). 

At At « 15-25 d and At > 35 d, the y- and J- 
band data are relatively sparse. We therefore calcu¬ 
late the bolometric luminosity for the entire light curve 
both including y and J {g'r'i'z'yJ) and excluding them 
{g'r'i'z') to determine the NIR to integrated flux ratio. 
We find that the fraction of flux at these wavelengths 
increases monotonically as a function of time, from 13% 
at At « 6 d to 23% at At « 29 d. For epochs when only 
g'r'i'z' observations were available, we add a fractional 
NIR contribution for the y- and J-band from our linear 
fit. At late times we adopt the last NIR ratio measure¬ 
ment at At « 29 d of ~ 23% instead of extrapolating our 
linear fit. This may underestimate the NIR contribution 
at late times. Our NIR contribution measurements are 
consistent with an a nalogous measurement for SN 2008D 
(Modjaz et al.|2009), tho ugh slightly smalle r than for SN 
2(J0ybb and SIN 20i0bh (Cano et al. 2011), which have 
maximum NIR contributions of 35-45% at At « 25 d. 

The resulting bolometric light curve is displayed in 
Figure Our associated uncertainty measurements in¬ 
corporate errors from the flux measurements, as well as 
bandpass uncertainties, but do not include errors intro- 
d uced from line a r interp olation. 


Lyman et al. 
collapse SIN bo. 


(2014b) have created a model for core- 
ometric corrections using two filters for 


and r' synthetic photometry. We ap ply the methodol¬ 
ogy described by Lyman et al. (2014b) and find that this 
model leads to excellent agreement with our brute-force 
g'r'i'z'yJ bolometric light curve. This confirms that we 
are not underestimating the ultraviolet and NIR contri¬ 
butions in our quasi-bolometric light curve. 

5. SPECTRAL ANALYSIS 
5.1. Isolating the SN Component 


nearby events. The corrections include ultraviolet and 
NIR contributions. Since SN 2013dx is at z = 0.145, 
we use our K-corrected spectra (see (5.1) to extract g' 


In a similar manner to that of § |4.1[ we wish to isolate 
the SN component from our spectroscopic observations 
of CRB 130702A. First, for absolute-flux calibration, we 
normalize our spectra to (uncorrected) broadband pho¬ 
tometry at the appropriate epoch. This accounts for slit 
losses caused by variable seeing. 

We next deredden our sp ectra of SN 2013dx in an anal¬ 
ogous manner to that of § |4.1| This includes contribu¬ 
tions from both the Milky Way [E{B — V) = 0.038 mag] 
and the host galaxy (Ay^host = 0.1 mag). 

To remove the afterglow contribution, we assume that 
the spectrum can be described at all times (and frequen¬ 
cies) as a power law of th e for m f^{t,a) oc with 

a = 1.25 and = 0.52 (§ 4.1). We normalize this func¬ 
tion to our (extinction-corrected) broadband photome¬ 
try, and subtract the appropriate power-law model at 
the epoch of each spectrum. 

Finally, we fit the LMI/DCT lat e-time photometry t o 
a variety of template galaxies from Kinney et al.| (|1996 ). 
Similar to D15, we find that the best-ht tempfate is a 
starburst galaxy, and we adopted this (appropriately nor¬ 
malized) as the host contribution to the spectra. We take 
this approach instead of using our final Keck/LRIS spec¬ 
trum because of the blue-shutter failure. 

After these corrections, only the SN component re¬ 
mains. Figure displays the resulting spectra of SN 
2013dx after smoothing. We exclude the hrst spectrum 
(At = 1.17d) because it is completely afterglow domi¬ 
nated. In addition, we manually excise nebular emission 
lines from the host galaxy of SN 2013dx. 

5.2. Comparison with Other Type Ic-BL SNe 

The early-time spectra of SN 2013dx are fairly fea¬ 
tureless, but after a week, broad (z; « 3 x 10^kms“^) 
features appear. Together with the lack of obvious H 
and He emission, this leads us to classify SN 2013dx as 
a broad-lined Type Ic SN (Ic-BL), as has been the case 
for essen tially all well-studied GR B-associated SNe thus 
far (e.g., Wo^sley fc Blooiii||2006 ). 

In Figure]^ we plot the spectrum of SN 2013dx ob¬ 
tained arouiid maximum light (At = 14.2 d) and at late 
times (At = 33.3 d) with pseudocontinuum removed, 
bandpass hltered, scaled, and binned to a common log¬ 
arithmic wavelength scale along with mean spectra of 
Type Ic-B L SNe both with and without GRBs from |Mod^ 
jaz et al. (2015). We note that the mean spectra include 
spectra of SJN 2013dx from D15, but this is one of many 
objects. The absorption features from mean Type Ic- 
BL SNe both with and without GRBs align well with 
SN 2013dx absorption features. This indicates that SN 
2013dx has similar photospheric velocities as other Type 
Ic-BL SNe. 

At maximum light SN 2013dx has a similar blueshifted, 
broad Si II A6355 line as both SNe Ic-BL with and with- 
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Fig. 7. — g'r'i'z'yj bolometric light-curve fit with Arnett-Valenti relation fit for 4 < t < 30d (solid black line). Bolometric light 
curves of other spectroscopically confirmed GRB-SNe are included for comparison. Note that these “bolometric” light curves are actually 
“quasi-bolometric” and cover different wavelength ranges. 


out GRBs. The blended Fe II absorption feature around 
a blueshifted wavelength of 4800 A is similar to that of 
SNe Ic-BL with GRBs but is weaker than that of SNe Ic- 
BL without GRBs. SN 2013dx has a stronger Ca II ab¬ 
sorption feature around 7900 A and weaker O I absorp¬ 
tion feature around 7200 A than most SNe Ic-BL both 
with and without GRBs. 

At later epochs, SN 2013dx does not deviate from the 
mean SN Ic-BL both with and without GRBs except be¬ 


yond 8200 A. However, the relative variation from the 
continuum seems weaker than the mean spectra of both 
SNe Ic-BL with and without GRBs. SN 2013dx has a 
similar Ga II absorption feature around 7900 A as SNe Ic- 
BL with GRBs but weaker than that of SNe Ic-BL with¬ 
out GRBs. 

To search for other similar objects in the literature, 
we use the cross-correlation tool SN Identification code 


(SNID; Blondin & Tonry 2007|. Several SNe Ic-BL 
that were not associated with GRBs, such as SN I997ef 


(Iwamoto et a l. 12000 ) and SN 20071 (Blondin et al. 12007 
Modjaz et al. 2014p , also provide good matches to SN 


2U13dx. In addition, DI5 hiRhl ight similarities to the en¬ 
ergetic SN 20I0ah (PTFlObzf; Corsi et al.|2011 Mazzali 
etahjlM^ . 


5.3. Photospheric Velocity Measurements 

In order to estimate the photospheric velocity of SN 
2013dx, we measure the velocity of the most prominent 
spectral feature, the Si II 6355 A absorption line. We 
employ a fitting code in IDL that removes a linear pseu- 
doc ontinuum and fits a Ga ussi an to the absorption l ine 
(see [Silverman et al. 2012 and Silverman et al. 2015 for 
a detailed description of tne code). i'able displays the 
inferred velocities for each spectrum. Our results are also 
consistent with those reported by D15. The first three 
spectra (At = 1.17, 3.25, and 6.22 d) are too noisy for 
reliable velocity measurements. 


We note that Parrent et al. (2015) suggests the ab¬ 
sorption feature at ^ 6100 A, normally identified as Si II 
6355 A, may be instead associated with Ha. Therefore, 
we compare mean SN Ic-B L photospheric veloc ities mea¬ 
sured using Fe II 5169A (Modjaz et al. 2015) at maxi¬ 


mum light. We find that our measurements are consis¬ 
tent with the measurements from the less-contaminated 
Fe II 5169 A. 


5.4. Line Identification 

We use SYN-I—I- ([Thomas et al.||2011|) to help identify 
the ions present in our spectra of SJN 2013dx. SYN-I--I- is 
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Fig. 8. — SN 2013dx spectrum with pseudocontinuum removed and bandpass filtered, binned to a common logarithmic wavelength 
scale (black). The mean (white) and standard deviation (red) of Type Ic-BL spectra from|Modjaz et al.|(|2015|l with (left) and without 
(right) GRB within 2 rest-frame days of the SN 2013dx spectrum, (top) Spectrum taken at At = 14.2 d, closest to maximum light in r' 
(At = 13.2d) and (At = 14.7d). (bottom) Spectrum taken at At = 33.3d. 


TABLE 4 

Velocity of Si II A6355 


At 

(d) 

Velocity 

(kms“^) 

9.3 

28,100 ± 1000 

11.3 

25,200 ± 500 

14.2 

21,300 ± 500 

31.3 

11,700 ± 500 

33.3 

10,800 ± 500 


Note. — In observer frame. The 
reported uncertainties come from fitting 
the Si II A6355 absorption feature with a 
single Gaussian function and do not in¬ 
clude errors from potential blending. 


derived from SYNOW (Fisher et al.|l99 

71), which uses the 

Sobolev approximation ( Sobolev||l96U 

Castor[ 

1970 

Jef- 


photospheric phase. SYN++ assumes that spectral lines 
are formed via resonance scattering above a sharp pho¬ 
tosphere. The location of the photosphere is expressed 
in velocity coordinates as Uph (in kms“^) and takes into 
account the homologous expansion of the ejecta. 

The optical depths for each species must also be input 
and line strengths are computed assuming Boltzmann 
excitation (i.e., local thermodynamic equilibrium, LTE) 
using a specified excitation temperature Tgxc (in K)- Non- 


LTE effects are partially taken in account by allowing 
different Texc values for each species, all of which can be 
different from the photospheric temperature Tphot- The 
latter is used only in computing the blackbody radiation 
emitted by the photosphere. All our SYN-I—I- fits are 
computed with all ions turned on simultaneously with a 
blackbody. 

We attempt to model the major spectral features of 
SN 2013dx at At = 9.3 and 33.3 d and look for evolution 
during the photospheric phase (see Figure [^. At 9.3 d 
after the burst, SYN-I--I- indicates a photos^eric veloc¬ 
ity of 30,000 kms“^, an outer velocity of the line-forming 
region of about 90,000 km s“^, and an estimated photo¬ 
spheric temperature of 16,000 K. The spectrum contains 
absorption from O I, Si II, and Fe II, and possible ab¬ 
sorption signatures of Fe III, Mg II, C II, Ca II, and 
Na I. We caution that because of the relatively uncer¬ 
tain line identifications, the derived velocities are robust 
indicators only for the lines of Fe II and Si II. 

By 33.3 d after the burst, according to our second 
SYN-I—I- fit, the SN ejecta have slowed down and cooled 
off significantly. This model indicates approximate val¬ 
ues for the photospheric velocity, outer velocity, and pho¬ 
tospheric temperature of 11,000kms“^, 60,000kms“^, 
and 9000 K, respectively. The majority of the absorp¬ 
tion in this spectrum is likely produced by Fe II and 
Ti II, though there is some evidence of Si II, Ca II, and 
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Fig. 9.— SYN++ fits to the 9.3 d (top) and 33.3 d (bottom) 
spectra of SN 2013dx. The spectrum of each individual ion is 
labeled. Their sum - plus a 16,000 K (top) and 9000 K (bottom) 
blackbody - is plotted in red on top of the actual observed spectra 
of SN 2013dx (binned to 6 A per pixel). 

possibly O I as well. While the fit to this spectrum at 
wavelengths below ^ 4700 A is not perfect, the broad 
peaks and troughs roughly match. This part of the spec¬ 
trum is notoriously difficult to model owing to hundreds 
of overlapping spectral features, mostly from iron-group 
elements. 


6. SUPERNOVA EXPLOSION PARAMETERS 
6 .1. Derived Parameters 

We model the basic explosion parameters of SN 2013dx 
by fitting its bolometr ic light curve wit h the Type I 
SN an alytical model of Arnett (1982) and Valenti et al. 
(2008 1 . This model assumes (1) homologous expansion of 
the ejecta, (2) spherical symmetry, (3) all ®®Ni is located 
at the center of explosion and no mixing, (4) radiation- 
pressure dominated ejecta, (5) the initial radius before 
explosion is small, ( 6 ) the diffusion approximation is ap¬ 
propriate for photons (i.e., the ejecta are in the photo- 
spheric phase), and (7) a single opacity over the duration 
of the explosion. 

The peak luminosity correlates with the mass of ®®Ni 


(Mni), while the light-curve shape is determined by the 
total ejecta mass (Mej) and the ejecta kinetic energy 
(Ak). We can break the degeneracy between Mej and 
Ak with photospheric-velocity measurements from our 
optical spectra. 

The timescale of the light curve is given by 


\j3cj J 


( 3 ) 


where /3 » 13.8 is an i ntegration constant. For a uniform 
density ( Arnett||l982[ p^ 




3 ^ejttph 
5 2 


( 4 ) 


We assume k = 0.07cm^g“^ to dire ctly compare with 
the literature for other GRB-SNe (e.g., Cano et al.j2011 ). 

We fit our quasi-bolometric light curve 575 ) with the 

Arnett-Valenti relation, 


Aph(t) =MNie 


px pX 

(cNi - eco) / A{z)dz + eco / B{z)dz 
^0 ^0 


( 5 ) 


with 


A{z) = 2ze-‘^^y+^\ 

B{z) = 

X = t/Tm, 

y = 7-m/(2TNi), and 

5 = 7'm(7'Co TNi)/(^TCoTNi)- 


( 6 ) 


The decay times of ®®Ni and ®®Co are tn, = 8.77d and 
Tco = 111.3 d, and the energies produced in one sec¬ 
ond by one gram of ®®Ni and ®®Co were taken as eNt = 
3.90 X 10^° ergs~^ g~^ and ec,n = 6.78 x lO^ergs"^ 


(Sutherland & Wheeler 1984 Cappellaro et al.| 1997). 


From our spectra and light curves, the SIN component 
was dominant starting at At « 4d (compare with SN 
2010bh, where sh ock breakout was prominent out to 7 d; 
Cano et al. 20111. The Arnett-Valenti relation assumes 
that the material is in the photospheric phase, which is no 
longer valid at At > 30 d. Therefore, our fit only includes 
4 < At < 30d. We find Mm = 0.37 ± 0.01 M© and 
Tm = 11.35±0.17d (statistical errors only). Usi ng Up h = 
21, 300kms~^ from our spectral fit near peak (§ |5.3[ ), we 
calculate Mej = 3.1 ± 0.1 Mq and Ak = (8.2 ± 0.43) x 

lO^i ej-g 


6 .2. Comparison with Other GRB-SNe 

We compare our bolometric light curve of SN 2013dx 
with that of other spectr oscopically confirme d CRB-SNe 
in Figure [7|(SN 1998 bw, Calama et al.|199^ SN 200 31w, 
[Mazzali eT al.||2006| SN 2U03dh, |Ueng et al.||2005[ SN 
2006aj, |Pian et al.|f2006[ SN 2009nz, [Olivares E. et al.| 

Note there is a typo incorrectly stating Ek ~ | — 2 "^ in 
the original text that was corrected W [Arnettl | |1996| |. This has 
been taken into account in Equations [3| and 1 his typo has been 
propagated throughout the literature. 
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TABLE 5 

Physical Parameters of GRB-SNe 


GRB-SN 

z 

-^• 7 , iso 

I'ph 

Mni 

Mej 

Ek 

Reference 



(erg) 

(km s~^) 

(Mo) 

(Mq) 

(1052 erg) 


GRB 980425/SN 1998bw 

0.0085 

(9.29 ± 0.35) X lO'^^ 

18,000 

0.42 ± 0.02 

6.80 ±0.57 

1.31 ± 0.10=" 

( 1 ), ( 2 ) 

GRB 030329/SN 2003dh 

0.1685 

1.33 X 10®2 

20,000 

0.54 ± 0.13 

5.06 ± 1.65 

1.21 ± 0.39=" 

( 1 ), ( 2 ) 

GRB 031203/SN 20031w 

0.105 

l-67lj( '(o ^ 

18,000 

0.57 ± 0.04 

8.22 ±0.76 

1.59 ±0.15=" 

( 1 ), ( 2 ) 

GRB 060218/SN 2006aj 

0.0335 

4.33l“ d ^ 

20,000 

0.21 ± 0.03 

2.58 ±0.55 

0.61 ± 0.14=" 

( 1 ), ( 2 ) 

GRB 091127/SN 2009nz 

0.49 

(4.3 ± 0.3) X 10®2 

17,000 

0.33 ±0.01 

4.69 ±0.13 

0.81 ± 0 . 02 =" 

(1), (3) 

GRB 100316D/SN 2010bh*' 

0.059 

> (3.9 ± 0.3) X lO'^® 

25,000 

0.12 ± 0.02 

2.47 ± 0.23 

0.92 ±0.08=" 

(1), (4) 

GRB 120422A/SN 2012bz'’ 

0.283 

4.5 X lO'^® 

20,500 

0.57 ± 0.07 

6.10 ±0.49 

1.53 ±0.13=" 

(1), (5) 

GRB 130427A/SN 2013cq 

0.3399 

(9.6 ± 0.04) X 1053 

32,000 

0.28 ± 0.02 

6.27 ± 0.69 

6.39 ± 0.70 

( 6 ) 

GRB 130702A/SN 2013dx 

0.145 

6.4l ( o X lO®® 

21,300 

0.37 ± 0.01 

3.1 ± 0.1 

0.82 ± 0.04 

(7) 

GRB 140606B/iPTF14bfu 

0.384 

(3.47 ± 0.02) X 1051 

19,820 

0.42 ±0.17 

4.8 ± 1.9 

1.1 ± 0.7=" 

(8) 


MgjuA 

Note. — ‘‘Ek originally calculated as Ek = - 2 “^’ scaled by a factor of 3/5 to directly compare with our values. All 

® 7 ,iso values calculated over 1 keV - lOMeV, except those indicated by superscript “b” which are calculated over 15-150 keV 
or superscript “c” wh ich are calculated over 2 0—1 400keV. (11 FCano l 1 201311, l2llKanek oet anil2007ll. (3l lTroiaet al.lll2012ll. (4) 
[Starling et al.| l |2011[ |, (5) [Zhang et al.| ( |2012| l, ( 6 ) [Xu et al.[}5jl3[ l, () this work, ( 8 j [Gano et al.[ ( [2015^ . 




Fig. 10.— (Isft) Comparison of SN explosion parameters from template GRB-SNe and SN 2013dx from Table|^ The size of the points 
indicates Upjj near the SN brightness peak. There does not appear to be any correlation between Mni and E-y^iso or Mni and Up^. (right) 
Comparison of SN explosion parameters using photometric data reported in the literature with our quasi-bolometric fitting procedure. This 
avoids using SN 1998bw as a template for other GRB-SN bolometric fits (see text for details). Blue points are fit to only optical data, red 
points are fit to optical and NIR data. Black line connects points that have both. 


2015[ SN 2010bh, [Olivares E. et al. 

[2012[ SN 2012bz: 

Meiandri et al.[[20I2 

and Schulze et a 

i.|2014| SN 2013cq: 

Meiandri et ai.||2014 

iPTFldbiu: Cano et al.||20I5 

). We 


nosity for SN 2012bz is assumed to be the same as that 
observed for SN 2010bh. Although these GRB-SN bolo¬ 
metric light curves cover different wavelength ranges, we 
can get a sense of the light-curve evolution. SN 2013dx 
most closely matches the light-curve shape of SN 2012bz, 
but with the caveat that SN 2013dx has a steeper rise 
than SN 2012bz. 


After 30 days, SN 2013dx appears to drop in luminosity 
rapidly. This is unlike the three bursts with late-time 
coverage - SN 20031w, SN 2003dh, and SN 1998bw - seen 
particularly well juxtaposed against SN 1998bw, which 
has extensive observations out to hundreds of days. This 
drop in luminosity is not from underestimating the NIR 
contribution at late times; we still observe a rapid drop 
in luminosity at 30 d if we continue the monotonically 


increasing NIR ratio function from § |4.3[ instead of adding 
a flat NIR contribution of 23% after 29 d. 

In order to compare the derived properties of SN 
2013dx with a broader sample of GRB-associated SNe, 
we use the derived values for Mn,, M r,, and Ek for all 
well-sampled events from Cano] (2013). The authors fit 
a template of SN 1998bw to deter mine the appropriate 
stretch and scale parameters (e.g., § 4.2). Using average 
s and k values for each burst, the authors then fit a scaled 
version of the UBVRIJH SN 1998bw light curve to the 
Arnett-Valenti model to derive the explosion parame¬ 
ters. This method has the benefit of (effectively) uniform 
wavelength coverage, even for events that were only ob¬ 
served in a few filters. However, the primary drawback is 
the assumption that all bolometric light curves are well 
fit by an appropriated scaled version of SN 1998bw. As 
evidenced by Figure this assumption breaks down at 
the very least for the redder filters for SN 2013dx (see also 
Lyman et al.||2014a|). Nonetheless, the explosion param- 
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Fi g. 11.— There is a clear correlation between explosion energy, Ek. and Mn-, for SNe lb, Ic . Ic-BL (|Cano|2013l and references therein), 
IIP ||Hamuy|2003| and references therein), and GRB-SNe from TableSee |Filippenko| l|1997|l for a review ot blN classification. 


eters for all well-studied GRB-associated SNe derived in 
this manneip^ are presented in Table 

To avoid any biases introduced by GRB-SNe that are 
not well matched with SN 1998bw, we also create our 
own quasi-bolometric light curves by compiling photom¬ 
etry from the literature (Table 1^. We fit these li ght 
curves with the Arnett-Valenti relation described in § |6.1| 
and break our results into those events with only optical 
datasets and those with optical and NIR photometry in 
Figure [l0| (right panel). This provides an estimate of the 
fundamental explosion parameters that is independent of 
any assumed similarity with SN 1998bw. 

We plot the inferred Mn; as a function of both the 
isotropic prompt gamma-ray energy release (Figure 


including other cor e-collapse ev ents . As shown by pre- 


and the derived SN kinetic energy (Ak; in Figure 11 


101 


6.3. Caveats and Conclusions 

From Figure it is clear that there is no correla¬ 
tion between Mni and the prompt energy release us¬ 
ing either method for calculating SN explosion param¬ 
eters of other GRB-SNe. For example, SN 2013dx has a 
comparable mass of synthesized ®®Ni as the subluminous 
GRB 980425 / SN 1998bw and the extremely luminous 
GRB13042 7A / SN 2013cq (| Levan et al.||2014t |Melandri 
et al.||2014|) . Even if we were to apply a beaming correc¬ 
tion, GRB 130427A would still have E~, several orders o f 
magnitude larger than GRB 980425 ( Perley et al.||2014 ), 
but comparable Mni- Similarly, there is no clear correla- 
tion between Mni and photospheric velocity at peak. Nu¬ 
merical simulations unambiguously predict that the mass 
of synthesized ®®Ni should be c orrelated with the degree 
of asymmetry in the explosion ( Gonzalez-Gasanova et al. 


2014 Umeda & Nomoto 2008); to the extent that our 
models faithfully reproduce the relevant (global) SN ex¬ 
plosion parameters, this result is clearly not borne out 
by the data. 

On the other hand, the explosion energy of the SN 
ejecta is clearly correlated with Mni, particularly when 

We note that |Can^ l |2013| l and |Cano et al.| | |2015[ l assume that 

E'k = —scaled their reported E^i values by a factor 
of 3/5. 


vious authors (e.g., Cano |2013 and Lyman et al. 2014a 
for recent compilations), GRB-associated SINe on aver¬ 
age have a higher mass of synthesized ®®Ni and larger 
kinetic energies than any other class of c ore-collapse SN e 
(except perhaps the superluminous SNe; Gal-Yam|2012 ). 
That said, the SN explosion energies are typically nar- 
rowly clustered and do not appear to significantly exceed 
1052 erg, consistent with (perhaps even indicative of) a 
magnetar origin for these events (|Mazzali et al. 2014). 
SN 2013cq (associated with GRB 130427A) appears to 
be a significant outlier in terms of its inferred Ak, how¬ 
ever, which remains to be fully explained. 

D15 created a bolometric light curve over the range 
3000-10,000 A extrapolated from u' and i'. We exam¬ 
ined the bolometric light curve from D15 and found 
that our peak luminosity is consistent with theirs when 
accounting for our additional NIR coverage. D15 re¬ 
port Mni « 0.2 Mq, a factor of two lower than our 
quoted value. On the other hand, they derive a total 
ejecta mass (Mgj « 7 ± 2 M 0 ) and a SN kinetic energy 
(Ek ~ (3.5 ± 1.0) X 10^^ erg) approximately a factor of 
two and four (respectively) larger than those presented 
here. D15 scale nu merical simulations of the similarly 
shaped SN 2003dh (|Mazzali et al.||2006|) to estimate SN 
explosion parameters, as opposed to using an analytical 
model (e.g., Arnett-Valenti) as adopted in this paper. 

Since Mni is closely related to the peak luminosity, we 
believe that our reported Mni estimate is more accurate 
than the value D15 report. When NIR contributions are 
included, SN 2013dx has a similar peak luminosity as SN 
1998bw. We therefore expect SN 2013dx to have a similar 
56Ni 

mass as SN 1998bw. Both numerical simulations 
and ana lytical models produce Mm » 0.4 Mq for SN 
1998bw ( Mazzali et al.|[2006 Gano 2013). 

The discrepancy m ejecta mass and kinetic ene rgy is 
caused rnainly by different opacity assumptions. |Maz-| 


zali et al. (|2006|) assumed an opacity, k = 0.5 Ye cm^ g“^, 
where Ye is the number of electrons per baryon. We 
assume the authors used Ye = 0.46 for iron and re¬ 
calculate our Arnett-Valenti fit. With an opacity of 
K = 0.02cm^g“^, we report Mej = 9.2 ± 0.2M© and 
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TABLE 6 

Photometry References for SN 1998bw-Independent Bolometric Light Curves 


GRB-SN 

Filter Coverage 

References 

GRB 980425/SN 1998bw 

UBVRIJHK 

Clocchiatti et al. I (1201 Ik 

GRB 030329/SN 2003dh 

UBVR 

jjeng et ai.| |j:^uuo|| 

GRB 031203/SN 20031-1^“^ 

— 


GRB 060218/SN 2006aj 

GRB 091127/SN 2009nz=" 

BVRIJHK 

Ferrero et al.| (|2006f ; |Kocevski et al.|l|2007f 

GRB 100316D/SN 2010bh'" 

grizJH 

Olivares E. et al.l d2012 i 

GRB 120422A/SN 2012bz 

grizj 

bcnulze et al.l (1261411 

GRB 130427A/SN 2013cq‘= 

BVRI 

Melandn et al.|lj2Ui4|l 

GKB 13U7U2A/yJN 2U13dx 

grizyj 

i nis work 

GRB 140606B/iPTF14bfu 

griz 

Cano et al.|(2015|l 


Note. — ^Simultaneous epochs near peak we re limited and were insu fficient for our analysi s 


Cano et al. 


120111 and 


Bufano et al. 


(2012 k 


tools. '^Used lower host-galaxy extinction from 
^Poor fit using our analysis tools; we fit to the bolometric curve provided by Hi. Fian (private 
communications). 


E’k = (2.5 ± 0.1) X 10®^ erg. We also note that from 
our fit of SN 1998bw in §4.2| we can see that the light- 
curve evolution of SN 1998bw does not match that of SN 
2013dx well (see Figure]^; hence, the different values of 
Mej and FIk are not nnexpected. 

We note that opac ity greatly affecte and i?K, but 
does not affect Mni. [Wheeler et al. (2015^ found that a 
conflict exists when comparing properties determined by 
fitting the peak (using Arnett-Valenti methods) to those 
deter mined from fitting the late-tim e tail (using methods 
from Clocchiatti & Wheeler 19971. This conflict is par¬ 
tially resblved~5j7usmgTmean opacity d etermined from 
both p eak and late-time tail parameters. [Wheeler et al. 


(2015) find that in general, the opacity is often overesti¬ 
mated in the literature and the typical mean opacity is 
K « 0.01 cm^ g~^. In view of this potential discrepancy, 
we report values for k = 0.01 cm^ g“^ (Mej « 21 Mq, 
Ek « 6 X 10®^ erg) and k = 0.1 cm^ g“^ (Mej « 2M0, 
Ek « 6 X 10®^ erg) to draw attention to the spread in 
Mej and Ak from this variable. 


driven SNe (e.g 

., Gonzalez-Gasanova et al. 

2014 

Maz- 

zali et al.[[2013( 

Umeda & JNomoto||20U8p imply that the 


The derived ejecta mass may therefore be biased by line- 
of-sight effects, and not representative of the total mass 
ejected in the explosion. 


7. DISCUSSION AND SUMMARY 

We present extensive optical and NIR photometry 
of GRB 130702A/SN 2013dx spanning 1-63 d after the 
gamma-ray trigger, and optical spectra covering 1-33 d 
after the trigger. At z = 0.145, GRB 130702A/SN 
2013dx is sufficiently close to clearly detect and model 
the underlying SN component that emerged a week after 
the burst. 

We isolate the SN component and present multi¬ 
band light curves, a quasi-bolometric {g'r'i'z'yJ) light 
curve, and spectra of SN 2013dx. Detection of the 
broad Si II A6355 absorption line at velocities approach¬ 
ing 3 X lO'^kms”^, combined with the absence of H 
and He features, indicates that SN 2013dx is a broad- 
lined SN Ic. We estimate the SN explosion parame¬ 
ters using the Arnett-Valenti analytical relation and in¬ 
fer Mm = 0.37 ± 0.01 Mq, Mej = 3.1 ± 0.1 M©, and 
Ak = (8.2 ± 0.4) X 10^1 erg. 


Our analysis allows us to compare SN 2013dx with 
other GRB-SNe, as well as other core-collapse SNe (those 
of identical spectral type and not). This is of particular 
interest because GRB 130702A is of intermediate E^^iso, 
between low-luminosity and cosmological GRBs. There 
seems to be no clear relation between Mni, Mej, or Ek 
with GRB isotropic energy (Figure 10), even when con¬ 
sidering beaming corrections. The SlN appears to not be 
imprinted with any information about the formation of 
the relativistic jet aside from the high photospheric ve¬ 
locity and lack of H and He that allows us to classify 
all GRB-SNe as Type Ic-BL. This is somewhat puzzling, 
given the predictions of a correlation between the degree 
of asymmetry and mass of synthesized ®®Ni for jet-driven 
explosions. We caution that the Arnett-Valenti relations 
we use to derive Mni assume spherical symmetry and 
this assumption may account for some but not all of the 
scatter in Figure [^ On the other hand, our observa¬ 
tions do provide support for predictions that Mni should 
be strongly correlated with the kinetic energy of the SN 
itself. 

Spectroscopically, SN 2013dx resembles both other 
GRB-SNe like SN 2006aj and SN 1998bw, as well as non- 
GRB SNe Ic-BL such as SN 1997ef, SN 20071, and SN 
2010ah. In terms of light curves, SN 2013dx most closely 
matches the evolution of SN 2012bz, associated with an 
intermediate GRB, but has a similar peak luminosity as 
SN 1998bw, associated with a low-luminosity GRB. Di¬ 
rect comparison of the light-curve evolution between SN 
2013dx and SN 1998bw indicates that SN 2013dx has a 
quicker rise time than SN 1998bw. The faster rise time 
may suggest that SN 2013dx has a steeper distribution of 
®®Ni in the o uter layers of the star (i.e., less mixing) t han 
SN 1998bw ([Piro & Nakar[[2013 fDessart et al.[[2012[). 

Late-time observations several months after the burst 


1 Wheeler et al. 

2015) and spectroscopy (Mazzali et al. 

»'Maeda et al.[2002 

Milisavljevic et al.|2015 

). We do 


conduct these effective asymmetry tests, but we strongly 
encourage late-time follow-up data for GRB-SNe when 
possible. 

Finally, we suggest two potential avenues for future 
study, especially with respect to GRB 130702A / SN 
2013d x. Detailed numeric al modeling of the SN ejecta 
(e.g., Mazzali et al. 2006 for SN 1998bw), specifically 
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tailored to the light curves and spectra of SN 2013dx 
(instead of simply scaling results from previous simu¬ 
lations), should help to improve the accuracy of esti¬ 
mates of the fundamental SN explosions parameters. In 
addition, a broadband study of the afterglow emission, 
in particular incorporating radio wavelengths, would en¬ 
able a much-improved estimate of the properties of the 
fastest-moving ejecta. This would greatly assist in plac¬ 
ing GRB 130702A in the context of other relativistic ex¬ 
plosions, specifically how the explosion ener gy is par¬ 


titioned wi th respect to ejecta velocity (e.g., Margutti 

eFaTjlMIl ). 
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TABLE 7 
Imaging log 


Filter 

Epoch 

(days) 

Telescope 

Exp. 

(s) 

AB mag 

u' 

330.49 

Keck/LRIS 

203 

24.50 ± 0.14 

g' 

1.17 

P60 

120 

18.80 ± 0.04 

g' 

1.26 

P60 

120 

18.86 ± 0.04 

g' 

2.21 

P60 

270 

19.52 ± 0.04 

g' 

3.23 

P60 

450 

20.02 ± 0.05 

g' 

4.19 

P60 

1260 

20.22 ± 0.04 

g' 

5.27 

P60 

540 

20.35 ± 0.04 

g' 

6.28 

P60 

540 

20.31 ± 0.04 

g' 

7.20 

P60 

720 

20.32 ± 0.04 

g' 

9.87 

Liverpool 

130 

20.35 ± 0.33 

g' 

10.87 

Liverpool 

130 

20.39 ± 0.35 

g' 

11.26 

P60 

720 

20.22 ± 0.05 

g' 

11.87 

Liverpool 

130 

20.06 ± 0.21 

g' 

12.18 

P60 

660 

20.23 ± 0.03 

g' 

12.86 

Liverpool 

130 

20.00 ± 0.23 

g' 

13.86 

Liverpool 

130 

20.39 ± 0.34 

g' 

14.18 

P60 

720 

20.31 ± 0.06 

g' 

15.18 

P60 

720 

20.43 ± 0.07 

g' 

16.18 

P60 

720 

20.36 ± 0.13 

g' 

17.87 

Liverpool 

150 

20.32 ± 0.13 

g' 

19.87 

Liverpool 

150 

20.78 ± 0.19 

g' 

21.87 

Liverpool 

150 

20.98 ± 0.15 

g' 

23.87 

Liverpool 

150 

21.25 ± 0.30 

g' 

27.92 

Liverpool 

150 

21.53 ± 0.07 

g' 

29.21 

P60 

540 

21.62 ± 0.12 

g' 

30.19 

P60 

1080 

21.72 ± 0.09 

g' 

34.17 

P60 

1080 

21.84 ± 0.08 

g' 

37.89 

Liverpool 

450 

22.06 ± 0.08 

g' 

39.88 

Liverpool 

450 

22.11 ± 0.09 

g' 

45.88 

Liverpool 

450 

22.64 ± 0.34 

g' 

53.88 

Liverpool 

150 

22.42 ± 0.07 

g' 

58.88 

Liverpool 

150 

22.56 ± 0.06 

g' 

62.87 

Liverpool 

150 

22.60 ± 0.09 

g' 

330.48 

Keck/LRIS 

203 

23.51 ± 0.03 

g' 

632.37 

DOT 

1200 

23.66 ± 0.04 

r' 

0.18 

P48 

60 

17.38 ± 0.04 

r' 

0.21 

P48 

60 

17.52 ± 0.04 

r' 

1.26 

P60 

360 

18.66 ± 0.05 

r' 

2.16 

RATIR 

7600 

19.23 ± 0.06 

r' 

2.21 

P60 

420 

19.32 ± 0.04 

r' 

3.21 

RATIR 

5040 

19.67 ± 0.04 

r' 

3.25 

P60 

120 

19.70 ± 0.08 

r' 

4.18 

P60 

1320 

19.97 ± 0.04 

r' 

4.21 

RATIR 

5360 

19.88 ± 0.05 

r' 

5.21 

RATIR 

4960 

19.95 ± 0.05 

r' 

5.27 

P60 

540 

20.00 ± 0.04 

r' 

6.20 

RATIR 

4880 

19.94 ± 0.05 

r' 

6.27 

P60 

540 

20.04 ± 0.04 

r' 

7.19 

RATIR 

2640 

19.90 ± 0.05 

r' 

7.19 

P60 

720 

20.06 ± 0.04 

r' 

8.22 

RATIR 

3920 

19.82 ± 0.05 

r' 

9.87 

Liverpool 

130 

19.86 ± 0.07 

r' 

10.87 

Liverpool 

130 

19.69 ± 0.07 

r' 

11.25 

P60 

660 

19.82 ± 0.02 

r' 

11.87 

Liverpool 

130 

19.73 ± 0.07 

r' 

12.18 

P60 

660 

19.78 ± 0.04 

r' 

12.22 

RATIR 

2880 

19.70 ± 0.04 

r' 

12.87 

Liverpool 

130 

19.74 ± 0.08 

r' 

13.16 

RATIR 

1280 

19.68 ± 0.05 

r' 

13.87 

Liverpool 

130 

19.68 ± 0.10 

r' 

14.18 

P60 

660 

19.70 ± 0.04 

r' 

14.18 

RATIR 

3760 

19.68 ± 0.07 

r' 

15.16 

RATIR 

5280 

19.70 ± 0.05 

r' 

15.17 

P60 

660 

19.64 ± 0.04 

r' 

15.87 

Liverpool 

150 

19.80 ± 0.07 

r' 

16.17 

P60 

660 

19.76 ± 0.08 

r' 

17.87 

Liverpool 

150 

19.83 ± 0.03 

r' 

19.87 

Liverpool 

150 

19.92 ± 0.04 

r' 

21.88 

Liverpool 

150 

20.00 ± 0.03 

r' 

23.87 

Liverpool 

150 

20.20 ± 0.02 

r' 

25.17 

RATIR 

4560 

20.16 ± 0.07 

r' 

26.15 

RATIR 

2160 

20.22 ± 0.05 
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TABLE 7 — Continued 


Filter 

Epoch 

(days) 

Telescope 

Exp. 

(s) 

AB mag 

r' 

27.15 

RATIR 

5200 

20.31 ± 0.05 

r' 

27.91 

Liverpool 

150 

20.42 ± 0.03 

r' 

28.16 

RATIR 

5120 

20.38 ± 0.05 

r' 

29.15 

RATIR 

5200 

20.52 ± 0.06 

r' 

29.19 

P60 

540 

20.49 ± 0.05 

r' 

30.18 

P60 

540 

20.54 ± 0.06 

r' 

30.19 

RATIR 

3120 

20.63 ± 0.06 

r' 

31.15 

RATIR 

3920 

20.62 ± 0.15 

r' 

32.15 

RATIR 

4960 

20.73 ± 0.06 

r' 

34.16 

P60 

540 

20.87 ± 0.07 

r' 

35.17 

RATIR 

3680 

20.96 ± 0.06 

r' 

37.88 

Liverpool 

150 

21.22 ± 0.04 

r' 

39.88 

Liverpool 

150 

21.27 ± 0.06 

r' 

41.16 

RATIR 

3840 

21.25 ± 0.08 

r' 

45.15 

RATIR 

2080 

21.33 ± 0.12 

r' 

45.87 

Liverpool 

150 

21.46 ± 0.14 

r' 

46.15 

RATIR 

3120 

21.36 ± 0.11 

r' 

47.88 

Liverpool 

150 

21.42 ± 0.15 

r' 

53.87 

Liverpool 

150 

21.75 ± 0.06 

r' 

58.87 

Liverpool 

150 

21.81 ± 0.06 

r' 

62.86 

Liverpool 

150 

21.80 ± 0.10 

R 

330.48 

Keck/LRIS 

200 

23.26 ± 0.06 

r' 

632.39 

DCT 

1200 

23.15 ± 0.04 

%' 

1.17 

P60 

120 

18.42 ± 0.04 

%' 

1.26 

P60 

120 

18.56 ± 0.06 

%' 

2.16 

RATIR 

7600 

19.09 ± 0.07 

%' 

2.21 

P60 

420 

19.12 ± 0.04 

%' 

3.21 

RATIR 

5040 

19.59 ± 0.06 

%' 

3.24 

P60 

450 

19.64 ± 0.05 

%' 

4.17 

P60 

1320 

19.97 ± 0.06 

%' 

4.21 

RATIR 

5360 

19.93 ± 0.06 

%' 

5.21 

RATIR 

4960 

20.05 ± 0.12 

%' 

5.26 

P60 

540 

20.13 ± 0.05 

%' 

6.20 

RATIR 

4880 

20.13 ± 0.09 

%' 

6.27 

P60 

660 

20.09 ± 0.05 

%' 

7.18 

P60 

660 

20.17 ± 0.05 

%' 

7.19 

RATIR 

3280 

20.07 ± 0.17 

%' 

8.22 

RATIR 

3920 

20.10 ± 0.07 

%' 

9.87 

Liverpool 

130 

19.93 ± 0.13 

%' 

10.87 

Liverpool 

130 

19.86 ± 0.16 

%' 

11.17 

P60 

180 

19.88 ± 0.11 

%' 

11.87 

Liverpool 

130 

19.87 ± 0.15 

%' 

12.19 

P60 

660 

19.94 ± 0.04 

%' 

12.22 

RATIR 

2880 

19.91 ± 0.05 

%' 

12.87 

Liverpool 

130 

19.94 ± 0.14 

%' 

13.16 

RATIR 

1280 

19.86 ± 0.05 

%' 

13.87 

Liverpool 

130 

20.01 ± 0.21 

%' 

14.18 

RATIR 

3760 

19.84 ± 0.05 

%' 

14.19 

P60 

660 

19.86 ± 0.06 

i' 

15.16 

RATIR 

5280 

19.83 ± 0.05 

i' 

15.19 

P60 

660 

19.77 ± 0.06 

i' 

15.87 

Liverpool 

150 

19.91 ± 0.10 

i' 

16.19 

P60 

660 

20.02 ± 0.10 

i' 

17.87 

Liverpool 

150 

19.94 ± 0.07 

i' 

19.87 

Liverpool 

150 

19.88 ± 0.05 

i' 

21.87 

Liverpool 

150 

19.97 ± 0.04 

i' 

23.87 

Liverpool 

150 

20.09 ± 0.05 

i' 

25.17 

RATIR 

4560 

20.04 ± 0.06 

i' 

26.15 

RATIR 

2160 

20.09 ± 0.05 

i' 

27.15 

RATIR 

5200 

20.15 ± 0.06 

i' 

27.91 

Liverpool 

150 

20.18 ± 0.02 

i' 

28.16 

RATIR 

5120 

20.21 ± 0.05 

i' 

29.15 

RATIR 

5200 

20.30 ± 0.07 

i' 

29.22 

P60 

540 

20.31 ± 0.06 

i' 

30.19 

RATIR 

3120 

20.41 ± 0.07 

i' 

30.20 

P60 

540 

20.44 ± 0.08 

i' 

31.15 

RATIR 

4000 

20.43 ± 0.14 

i' 

32.15 

RATIR 

5440 

20.52 ± 0.07 

i' 

34.18 

P60 

540 

20.67 ± 0.07 

i' 

35.17 

RATIR 

3680 

20.68 ± 0.07 

i' 

37.88 

Liverpool 

150 

21.01 ± 0.05 

i' 

39.87 

Liverpool 

150 

21.03 ± 0.08 

i' 

41.16 

RATIR 

3840 

20.88 ± 0.08 
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TABLE 7 — Continued 


Filter 

Epoch 

(days) 

Telescope 

Exp. 

(s) 

AB mag 

i' 

45.15 

RATIR 

2080 

21.05 ± 0.11 

i' 

45.87 

Liverpool 

150 

21.16 ± 0.10 


46.15 

RATIR 

3120 

21.12 ± 0.10 

i' 

47.88 

Liverpool 

150 

21.14 ± 0.13 


53.86 

Liverpool 

150 

21.56 ± 0.09 

i' 

58.86 

Liverpool 

150 

21.65 ± 0.09 


62.86 

Liverpool 

150 

21.53 ± 0.12 

i' 

632.39 

DOT 

1200 

23.03 ± 0.06 

z' 

2.16 

RATIR 

5700 

18.91 ± 0.09 

z' 

2.21 

P60 

870 

19.02 ± 0.08 

z' 

3.21 

RATIR 

3780 

19.36 ± 0.07 

z' 

3.24 

P60 

840 

19.43 ± 0.09 

z' 

4.19 

P60 

1320 

19.53 ± 0.10 

z' 

4.21 

RATIR 

4020 

19.72 ± 0.12 

z' 

5.20 

P60 

1020 

19.84 ± 0.12 

z' 

5.21 

RATIR 

3720 

19.84 ± 0.08 

z' 

6.20 

RATIR 

3660 

19.94 ± 0.11 

z' 

6.26 

P60 

540 

19.95 ± 0.12 

z' 

7.18 

P60 

540 

20.04 ± 0.11 

z' 

7.19 

RATIR 

2820 

19.71 ± 0.35 

z' 

8.22 

RATIR 

2940 

19.98 ± 0.13 

z' 

9.87 

Liverpool 

150 

20.11 ± 0.12 

z' 

10.87 

Liverpool 

150 

19.99 ± 0.12 

z' 

11.87 

Liverpool 

150 

20.13 ± 0.11 

z' 

12.22 

RATIR 

2160 

19.94 ± 0.13 

z' 

12.87 

Liverpool 

150 

20.02 ± 0.11 

z' 

13.16 

RATIR 

960 

20.02 ± 0.12 

z' 

13.87 

Liverpool 

150 

19.82 ± 0.13 

z' 

14.18 

RATIR 

2820 

19.87 ± 0.12 

z' 

15.16 

RATIR 

3960 

20.01 ± 0.10 

z' 

15.87 

Liverpool 

150 

20.21 ± 0.20 

z' 

17.88 

Liverpool 

150 

20.10 ± 0.08 

z' 

19.88 

Liverpool 

150 

20.08 ± 0.05 

z' 

21.88 

Liverpool 

150 

20.19 ± 0.06 

z' 

23.88 

Liverpool 

150 

20.27 ± 0.04 

z' 

25.17 

RATIR 

3300 

20.46 ± 0.14 

z' 

26.15 

RATIR 

1560 

20.06 ± 0.13 

z' 

27.15 

RATIR 

3900 

20.30 ± 0.09 

z' 

27.90 

Liverpool 

150 

20.33 ± 0.07 

z' 

28.16 

RATIR 

3780 

20.23 ± 0.09 

z' 

29.15 

RATIR 

3900 

20.49 ± 0.10 

z' 

30.19 

RATIR 

2340 

20.44 ± 0.15 

z' 

31.15 

RATIR 

2940 

20.45 ± 0.11 

z' 

32.15 

RATIR 

4080 

20.81 ± 0.15 

z' 

34.15 

RATIR 

3660 

20.87 ± 0.12 

z' 

35.17 

RATIR 

2820 

20.96 ± 0.14 

z' 

37.87 

Liverpool 

150 

21.15 ± 0.15 

z' 

39.87 

Liverpool 

150 

21.33 ± 0.20 

z' 

45.86 

Liverpool 

150 

21.41 ± 0.30 

z' 

53.85 

Liverpool 

150 

21.27 ± 0.13 

z' 

58.85 

Liverpool 

150 

21.60 ± 0.20 

z' 

62.85 

Liverpool 

150 

21.85 ± 0.29 

z' 

632.41 

DCT 

1200 

23.02 ± 0.12 

y 

2.16 

RATIR 

5700 

18.78 ± 0.08 

y 

3.21 

RATIR 

3780 

19.23 ± 0.08 

y 

4.21 

RATIR 

4020 

19.53 ± 0.09 

y 

5.21 

RATIR 

3720 

19.70 ± 0.09 

y 

6.20 

RATIR 

3660 

19.79 ± 0.07 

y 

12.22 

RATIR 

2160 

19.49 ± 0.10 

y 

13.16 

RATIR 

960 

19.70 ± 0.13 

y 

14.18 

RATIR 

2820 

19.75 ± 0.09 

y 

15.16 

RATIR 

3960 

19.73 ± 0.11 

y 

22.16 

RATIR 

2340 

19.85 ± 0.13 

y 

25.17 

RATIR 

3300 

19.69 ± 0.14 

y 

26.15 

RATIR 

1560 

20.00 ± 0.13 

y 

27.15 

RATIR 

3900 

20.05 ± 0.09 

y 

28.16 

RATIR 

3780 

20.01 ± 0.09 

y 

29.15 

RATIR 

3900 

20.03 ± 0.10 

y 

30.19 

RATIR 

2340 

20.04 ± 0.11 

y 

31.15 

RATIR 

2940 

20.28 ± 0.12 

y 

32.15 

RATIR 

4080 

20.00 ± 0.11 

y 

34.15 

RATIR 

3660 

20.49 ± 0.14 

y 

35.17 

RATIR 

2820 

20.51 ± 0.14 

j 

2.16 

RATIR 

6380 

18.75 ± 0.09 
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TABLE 7 — Continued 


Filter 

Epoch 

(days) 

Telescope 

Exp. 

(s) 

AB mag 

J 

3.21 

RATIR 

4230 

19.22 ± 0.07 

J 

4.21 

RATIR 

4500 

19.65 ± 0.09 

J 

5.21 

RATIR 

4160 

19.64 ± 0.10 

J 

6.20 

RATIR 

4090 

19.85 ± 0.10 

J 

14.18 

RATIR 

3150 

20.07 ± 0.12 

J 

15.16 

RATIR 

4430 

19.86 ± 0.11 

J 

27.15 

RATIR 

4360 

19.92 ± 0.13 

J 

28.16 

RATIR 

4230 

20.05 ± 0.12 

J 

29.15 

RATIR 

4360 

20.23 ± 0.13 

J 

349.45 

Keck/MOSFIRE 

132 

23.18 ± 0.32 

H 

2.16 

RATIR 

6380 

18.60 ± 0.09 

H 

3.21 

RATIR 

4230 

18.96 ± 0.08 

H 

4.21 

RATIR 

4500 

19.37 ± 0.11 

H 

5.21 

RATIR 

4160 

19.83 ± 0.12 

H 

6.20 

RATIR 

4090 

19.89 ± 0.13 

Ks 

349.44 

Keck/MOSFIRE 

79 

> 22.29 


Note. — AB magnitudes, not corrected for Galactic extinction. 



